*Table C.1 on correlates of data missingness*

version 15.1

clear

use "worms_replication.dta"
set more off


***create interactions with each measure of wealth***
gen ConcTrade = Concordats*trade_fixed
gen ConcCalorie = Concordats*calories
gen Conc_pop_5 = Concordats*city_5_fixed
gen Conc_pop_10 = Concordats*city_10_fixed
gen Conc_pop_15 = Concordats*city_15_fixed

	
*Create half Century Counters*
gen half_century = 50*floor(year/50)
label var half_century "Half Century"
quietly forval i = 3/18{
	gen HC`i'0 = 0
	gen HC`i'5 = 0
	replace HC`i'0 = 1 if year >= `i'*100 -50 & year <`i'*100
	replace HC`i'5 = 1 if year >= `i'*100 & year <`i'*100 + 50
	}


	
/* Collapse data to diocese-bishop observations */
sort diocese year
drop if byearin != year 


***generate missingness measure***
gen missing = 0
replace missing = 1 if religious_bishop == .
label var missing "Bishop Alignment Unobserved"


***Table C.1, each regression is a column***
reghdfe  missing ConcTrade Concordats crusade1 crusade2 crusade3 crusade4 if  year<=1309,  a(iddiocese half_century) vce(cluster iddiocese)
reghdfe  missing ConcCalorie Concordats crusade1 crusade2 crusade3 crusade4 if  year<=1309,  a(iddiocese half_century) vce(cluster iddiocese)
reghdfe  missing Conc_pop_5 Concordats crusade1 crusade2 crusade3 crusade4 if  year<=1309,  a(iddiocese half_century) vce(cluster iddiocese)
reghdfe  missing Conc_pop_10 Concordats crusade1 crusade2 crusade3 crusade4 if  year<=1309,  a(iddiocese half_century) vce(cluster iddiocese)
reghdfe  missing Conc_pop_15 Concordats crusade1 crusade2 crusade3 crusade4 if  year<=1309,  a(iddiocese half_century) vce(cluster iddiocese)
